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Abstract 

For relativistic particles with spin 1/2, which are described by the Dirac equation, 
a semiclassical trace formula is introduced that incorporates expectation values of 
observables in eigenstates of the Dirac-Hamiltonian. Furthermore, the semiclassical 
limit of an average of expectation values is expressed in terms of a classical average 
of the corresponding classical observable. 
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1 Introduction 

The three decades that passed since the completion of Gutzwiller's pioneering work on the 
trace formula |l|] (see also 0) have witnessed a flourishing development of semiclassical 
methods in connection with trace formulae. Gutzwiller himself, e.g., applied his semiclassi- 
cal representation of the spectral density in terms of a sum over classical periodic orbits to 
the quantisation of the anisotropic Kepler problem [|1], |3j . Later, in the context of the then 
developing field of quantum chaos, it was realised that the Gutzwiller trace formula can 
quite generally be employed for a semiclassical quantisation of classically chaotic systems, 
where no Bohr-Sommerfeld (or Einstein-Brillouin-Keller) type quantisation is available. A 
number of semiclassical quantisation schemes were invented and successfully passed vari- 
ous tests, see 0, |J |] for collections of examples. On the other hand, Hannay and Ozorio 
de Almeida and Berry |7j realised that a semiclassical analysis of spectral correlations 
can be based on trace formulae. These techniques, which have meanwhile been refined in 
several aspects, still provide the strongest analytic support for the conjectured universality 
in energy level statistics ||. 

Shortly after Gutzwiller established his trace formula, and apparently independently 
thereof, mathematicians began to prove various versions of semiclassical trace formulae 
B 1^' 11' 12' HI- Moreover, generalisations of the trace formula as a semiclassical repre- 



sentation of the spectral density were developed, e.g., to include also expectation values of 
observables |14], |P2f . Thus, by now the Gutzwiller trace formula is a well developed and 
well established tool in semiclassical quantum mechanics. 

So far, trace formulae have been derived for quantum systems with only translational 
degrees of freedom. These possess a well defined classical limit and thus semiclassical 
methods can be applied in a straight forward manner. From a physical point of view, 
however, spin degrees of freedom are of particular importance. Quantum mechanically 
a spin 1/2 of, say, an electron is described in a two- (i.e. finite-) dimensional Hilbert 
space, lacking a direct classical counterpart, so that a semiclassical description of quantum 
systems with coupled translational and spin degrees of freedom is less obvious than in 
the previous cases. However, the semiclassical formalism that was primarily developed in 
the mathematical community (see, e.g., JTD], [TTj, [T2|]) can also be applied to derive a trace 
formula for Dirac- and Pauli-Hamiltonians JL5J, jl6| . For these, alike in Gutzwiller's original 
trace formula, the spectral density can be expressed in terms of a sum over classical periodic 
orbits (of the translational motion). Spin enters through weight factors that depend on a 
phase derived from a 'classical' spin precessing along the classical periodic orbits. 

It is the goal of this paper to develop a semiclassical theory for expectation values of 
observables in the case of relativistic particles with spin 1/2. To this end we first extend the 



trace formula [ I5[ |T6f for Dirac-Hamiltonians to also include expectation values and then 
discuss averages of expectation values in eigenstates of the Hamiltonian that correspond 
to a certain spectral interval. In this way average expectation values and fluctuations 
thereabout are expressed in terms of classical quantities. The article is organised as follows: 
In section 2 we recall spectral properties of Dirac-Hamiltonians and represent them as 
Weyl operators. Bounded observables are then introduced in section 3. Subsequently, 



in section 4 the semiclassical description of the time evolution operator in terms of a 
Van Vleck-Gutzwiller propagator for Dirac-Hamiltonians that was developed in [|15|, [16| is 
discussed. Section 5 is then devoted to the derivation of the trace formula for the spectral 
density weighted by expectation values. Finally, in section 6 we calculate averages of 
expectation values semiclassically. 

2 The Dirac-Hamiltonian 

In relativistic quantum mechanics a point particle of mass m and charge e that is exposed 
to external electromagnetic fields E and B is described by the Dirac equation. Here we 
fix an inertial frame and work with fixed coordinates (t, x) for Minkowski space-time. We 
moreover assume that in the given frame the fields are static so that the Dirac equation 
reads 

ith^(t,x) = H4>{t,x), (2.1) 



with the Dirac-Hamiltonian 



H = ca- (- V - -A(xj) + mc 2 (3 + eip{x) , (2.2) 

where <p, A are electromagnetic potentials so that E(x) = —V(p(x) and B(x) = V x A{x). 
The Dirac algebra a^a u + a u a^ = 25^, //, v = 0, . . . , 3, is realised by the hermitian 4x4 
matrices 

cr\ , n [t 2 



a= {<r O) and ao = (3= {o -t 2 ) ' (2 - 3) 

where o^, k = 1,2,3, are the Pauli matrices and 1„ denotes the n x n unit matrix. For 
more details see [|l7j . 

If the potentials </?, A^ are smooth functions, or at most have Coulomb singularities at 
some points Xj such that 

|e^(a;)|, \eA k (x)\ < ^ + b (2.4) 

holds for some constants a < y, b > and for all x E M?\{xj}, the Dirac-Hamiltonian 
( p.2|) is essentially self-adjoint on the domain C^(M 3 \{x :; }) Cg> C 4 in the Hilbert space TC = 
L 2 (IR 3 ) (8> C 4 of square- integrable Dirac spinors (see, e.g., ||17|| ). The quantum Hamiltonian 
is then given by the (unique) self-adjoint extension which we also denote by H. As an 
example consider an electron in the Coulomb potential of a nucleus with charge Ze, such 

T 2 2 

that eip(x) = — -r|-. Since the fine structure constant a = f- is approximately given by 
j|y, the bound ( |2.4j ) implies essential self-adjointness for Z < 68. In the case A = and 
<£> spherically symmetric, however, the bound on a can be improved by a factor of v3, 



leading to Z < 118. For larger Z one has to choose special self-adjoint extensions, see |]T7 
for further details. 



Dirac-Hamiltonians ( |2.2| ) possess spectra that are not bounded from below. In fact, if 
the electromagnetic fields vanish as \x\ — > oo, the essential spectrum of H is given by 

a ess (H) = (-00, -mc 2 } U [+mc 2 , +00) . (2.5) 

(This is also true under appropriate weaker conditions, see ||17| .) In what follows we will 
always consider situations where H is (essentially) self-adjoint and where (|2.5|) holds. Since 
below we are only interested in properties of the discrete spectrum of Dirac-Hamiltonians, 
we will therefore concentrate on the eigenvalues E n of H with \E n \ < mc 2 . More precisely, 
we choose some energy E in the gap of a ess (H) and introduce an interval 

I(E; ft) := [E - fuv, E + ftw] , uj > fixed, (2.6) 

that has no overlap with the essential spectrum of H. This only requires h to be sufficiently 
small. We are then interested in the eigenvalues E n e I{E; ft) and the expectation values 
{ipn, Bip n ) of observables B in eigenstates ip n related to these eigenvalues. 

For the purpose of semiclassical asymptotics it is advantageous to view the Hamiltonian, 
as well as all further observables, as Weyl operators. For simplicity we now assume that 
the potentials are smooth so that the Weyl representation 

{Hi/>)(x) = 777^ / / H(p,^±y)er^-y)^(y) dpdy (2.7) 

{i-Khy j R3 j R 3 v 2 / 

of H holds for ip e C^(R 3 ) ® C 4 . Here 

H(p, x) := col- (p A(x) ) + mc 2 j3 + e(p(x) 



f(e^{x) + mc 2 )± 2 (cp — eA(x)) ■ a 
\(cp — eA(x))-cr (e<f (x) — mc 2 ) I2 



(2.8) 



is the Weyl symbol of H , which for all (p, x) e M 3 x R 3 is a hermitian 4x4 matrix with 
the two twofold degenerate eigenvalues 

#*(?, x) = &f{x) ± v/(cp-eA(x)) 2 + m 2 c 4 , (2.9) 

such that H + — H > 2mc 2 > 0. These eigenvalues are the classical Hamiltonians for 
relativistic particles (without spin) with positive and negative kinetic energies, respec- 
tively, that are exposed to the electromagnetic fields generated by the potentials tp and 
A. The presence of these two eigenvalues is a remnant of the fact that the Dirac equation 
( p.l|) describes both particles and anti-particles, and the dimension of the corresponding 
eigenspaces follows from the spin 1/2 of the particles and anti-particles, respectively. 

One can now introduce an orthonormal basis of C 4 consisting of the eigenvectors 
e\{p, x), e 2 (p, x) of H (p, x) with eigenvalue H + (p, x) and the eigenvectors /i(p, x), f2(p, x) 



with eigenvalue H (p,x). This basis allows to construct the 4x2 matrices V + (p, x) and 
V-(p, x) with ei(p,x),e 2 (p,x) and fi(p,x), f 2 (p,x), respectively, as columns [15|, |16|| , 

V + (p,x)= l f(e(p,x)+m^)t 2 



V-(p,x) 



2s(p, x){e(p, x) + mc 2 ) \ (cp - eA(x)) ■ a- 
1 / {cp — eA{x)) ■ <j 



(2.10) 



v/2e(p,a;)(e(p,a;) + mc 2 ) \-(^P, x ) + ™c 2 )l 2/ 

where 

s(p, x) := y/(cp - eA(x)) 2 + m 2 c 4 . (2.11) 

The columns of V± being orthonormal implies that these matrices represent isometries from 
C 2 into C 4 so that V±V± = t 2 - Moreover, the projectors P^(p,x) onto the eigenspaces 
corresponding to the eigenvalues H ± (p, x) are given by P$ = V±V±. 

3 Observables 

The kind of observables we are going to consider are given by matrix valued Weyl operators 
as in ( p.7|) , i.e. 

{B^){x) = — Lj / / B(p,^±y-,h)e^< x -y^(y) dpdy , (3.1) 

{2ti n) 6 J R 3 J t s V 2 / 

that yield bounded and self-adjoint operators on Ti. Notice that here the symbol B(p, x; H) 
is allowed to depend on h. In order to ensure that the operator B defined by ( |3.1| ) is 
bounded one can impose the following condition on its symbol: it shall be smooth in the 
variables (p, x) and, moreover, possess an asymptotic expansion 

B(p 7 x;H)~J2fr JB AP,x) (3-2) 

in the sense that for all iV > and all multi-indices a, (3 G Nq 

N 

mdi (B(p,x;h) -YttfBifavJ) < C aP h N+1 . (3.3) 

V z — ' / 4x4 

j=0 

In particular, B and all of its derivatives shall be bounded. Here || . . . || 4x4 denotes an 
arbitrary matrix norm and the short-hand d% = d" 1 " a 2 2 a3 a3 \ ias been used. The leading 

0X-i 0Xi-) 0Xrt 



term Bq(p,x) in the semiclassical asymptotics ( |3.2|) is called principal symbol and will be 
viewed as the classical observable corresponding to the quantum observable B. For scalar 
symbols a proof for the fact that (|3.3|) implies the boundedness of B can, e.g., be found 
in |Tj|; this immediately carries over to the present situation, see also [I5|. The operator 



B defined by ( |3.1| ) for ip e C^°(R 3 ) (g> C 4 hence can now be extended to all of 7i and is 
self-adjoint as soon as its symbol B(p,x; h) is hermitian. In the following we will always 
assume this to be the case. 

Two particular classical observables are the projectors Pq (p, cc) onto the eigenspaces 
of the symbol ( |2.8| ) of the Dirac-Hamiltonian. If we assume that the potentials and all 
of their derivatives are bounded, which poses no severe restriction since we want the elec- 
tromagnetic fields to vanish as \x\ —* oo, an inspection of fl2.10|) reveals that P^(p,x) 



together with all derivatives are also bounded. A Weyl quantisation of P^(p,x) hence 
yields bounded operators Pq. These, however, are not yet projection operators on TC, 
but a semiclassical construction, see pO, [H], allows to introduce bounded and self-adjoint 



operators P ± with principal symbols P (p, x) that obey P + + P~ = 1%. Moreover, up to 
terms of order h°° these operators fulfill (-r) 2 = P , P + P~ = and commute with H. 
With respect to these semiclassical projectors one can now split an observable B of the 
type introduced above according to 

B = B d + B nd} (3.4) 

where 

B d := P + BP + + P-BP- and B nd := P + Bp- + P~BP + (3.5) 

denote the diagonal and the non-diagonal part, respectively. According to the product rule 
for Weyl operators (see, e.g., [|^, |i~8| ) B d and B nd are again bounded and self-adjoint Weyl 
operators with principal symbols 

B d ,o(p, x) = P + (p, x) B (p, x) P + (p, x) + Pq(j>, x) B (p, x) P (p, x) , 
Bndfi{p, x) = P + (p, x) B (p, x) P ~(p, x) + P (p, x) B (p, x) P+(p, x) . 

4 Semiclassical time evolution 

In this section we mainly review the derivation of the Van Vleck-Gutzwiller propagator 
for Dirac-Hamiltonians that was developed in [[T|, [RJ. For a detailed exposition see [I6|. 



The aim here is to obtain the semiclassically leading asymptotics of the Schwartz kernel 
K(x,y,t) for the unitary time evolution operator U{t) = e~n tH . Given an initial state 
■00 ^ W> the solution ip(t) = U(t)il) of the Dirac equation (|2.1| ) with initial condition 
i/j(0) = ipo reads in coordinate representation 

V>(t,a?)= / K(x,y,t)ip {y) dy , (4.1) 

so that the initial condition for the kernel is K(x, y, 0) = 5(x — y). Up to terms of order 
h°° this kernel can be approximated by a semiclassical Fourier integral operator, i.e. 

K(x, y, t) = j-L- [ \a+(x, y, t, e^^ 

+a^{x,y,t,€)e^-( x >y' t >®] d^ + 0{h°°) . 



The two additive contributions under the integral refer to the two eigenvalues H ± of the 
symbol H of the Dirac-Hamiltonian, which each serve as a classical Hamiltonian. The 
amplitude factors a h take values in the 4x4 matrices and are assumed to allow for the 
asymptotic expansions 

a±(x,y,t,$) = J2{j) "kfaV't'Z) > ( 4 - 3 ) 



fc>0 



with initial conditions 




at{x,v,0,£) = t»™ ' . ff ' (4.4) 



Together with the initial condition ± |t=o = $, ■ (x — y) for the phase function this ensures 
the desired behaviour of the representation Q4.2 ) for t — > 0. 



The strategy to obtain a semiclassical approximation for K(x } y, t) now is the same as 
in the WKB method: one inserts the ansatz Q4.2Q with the expansion ( |4.3|) into the Dirac 
equation, groups terms together with like powers of h, and equates these terms to zero. 
This yields a hierarchy of equations for the phase functions <jy^ and the coefficients a k that 
have to be solved order by order in k. The lowest order (k = 0) equation, together with 
the definition S ± (^,x,t) := ± (a3, y, t,£) + £ ■ y, relates the phase functions ^ to the 
solutions S ± of the associated Hamilton- Jacobi equations 

f)Q± 

H ± (V x S ± (£,x,t),x) + —(£,x,t) = with S ± (tx,0) = £-x . (4.5) 

Hence S ± generate canonical transformations (p ± ,x) h-> (^, z^) with p^ = V x 5 ,± (^, x, t) 
and z ± = V^5' ± (^, x, t) such that (^,2 ± ) and (p ± ,x) are starting and end points, re- 
spectively, of solutions (P ± (t'), X ± (t / )), < t' < t, of the classical equations of motion 
generated by H ± . 

The leading semiclassical asymptotics for the amplitudes then follows from both the 
equations in leading and next-to-lowest order. For sufficiently small t the result reads 



V±(V x S±(Z,x,t),x)d ± (V x S±(£,x,t),x)Vi(£,y) . 

The restriction to small t stems from the non-uniqueness of the solutions S^ 1 of the 
Hamilton- Jacobi equations (|4.5| ) for times beyond which caustics appear. One is, however, 
faced with the same problem already in the context of non-relativistic, spinless particles 
where the respective solution merely consists of the first factor on the right-hand side of 
( 3~E ). For this case Gutzwiller showed [E3| that one can pass beyond caustics with a correct 



choice of the phase in the first term of (OS) when sign changes occur under the square-root. 



The appropriate phase factor that Gutzwiller introduced essentially consists of the Morse 



index of the related classical trajectory. A mathematically rigorous justification of this 
construction can, e.g., be found in 



10 



|TTJ. Obviously, exactly the same procedure can be 
applied to extend (|4.6| ) to arbitrary, however finite, times. 

The factor d± appearing in ( |4.6|) is a solution of the spin transport equation 

d 
At 



-iM ± (P ± (t),X ± (t)) 



d ± (P ± (t),X ± (t))=0 



(4.7) 



with initial condition d±(P (0), X (0)) = I2, in which the time derivative is to be taken 
along the solutions (P (£), X (£)) of the classical equations of motion that are generated 
by the solutions S^ of the Hamilton- Jacobi equations ( }4.5| ). The 2x2 matrices 



M ± (p,x):=T 



ec 



2e(p,x) 



B(x)± 



e(p, x) + mc 



(e{x)x(p-- c A{x)) 



■ (j 



(4- 



are hermitian and traceless so that according to ( f4.7|) the spin transport matrices d± take 
values in SU(2). Since the solutions (P (t),X (£)) of Hamilton's equations of motion 
depend on the initial points (£, z^), below we also write d±{£, z ± 1 1) for d±(P ± (t), X ± (t)). 
In these variables the composition law of the spin transport matrices reads 



d ± (Z,z ± ,t + t')=d ± (P ± (t'),X ± (t'),t)d ± (ti,z ± ,t') 



(4.9) 



In order to calculate the leading semiclassical asymptotics of the time evolution kernel 
K(x, y, t) further one employs the method of stationary phase to the integral (|4.2|) . Due to 
the connection of the phases ± with the solutions S ± of the Hamilton- Jacobi equations, 
the stationary points £ st of ± are such that the canonical transformations generated by 
S ± (£ st , x, t) map (p ± ,x) to (£ si ,y). These are therefore the end and starting points, 
respectively, of the associated solutions (P ± (t / ), X ± (t')) of the equations of motion. Thus 
the stationary points £ st are in one-to-one correspondence with the classical trajectories 
7 1 * 1 connecting y to x in time t. The expression that finally results from the application of 
the method of stationary phase is to a large extent analogous to the Van Vleck-Gutzwiller 
propagator (see |25|) for non-relativistic, spinless particles. Differences arise due to (i) the 
presence of two kinds of classical dynamics (generated by the two Hamiltonians H^) instead 
of one, and (ii) the presence of the additional factors in the amplitude fl4.6p . The latter, 
however, need only be evaluated at the stationary points £ st . As a result one therefore 
obtains 



K(x,y,t) 



(2vrin) 3 / 2 



J2 V ±(t)d±Vl(0)D, 



± e 



k R ~±- 



■ 1 ?",± 



(1 + 0(h)) 



(4.10) 



Here the sum extends over all solutions 7 of the equations of motion from y to x in 
time t, and z/ 7 ± denotes their Morse indices. Moreover, the phase functions ^ evaluated 
at the stationary point $, st related to 7 ± turn out to be Hamilton's principal functions 



R±(x,y,t) = S±(ti st ,x,t)-t 



St 



y 



L> 7 ± 



and 



det 



d 2 R\ 
dxdy 



{x,y,t) 



1/2 



(4.11) 



We also introduced the abbreviation V±(tf) := V ± (P ± (t'),X ± (t')). 

Inspecting ( 4.10|) reveals a close similarity to the Van Vleck-Gutzwiller propagator. 



The summation goes over the classical trajectories of relativistic particles, such that here 
only the translational dynamics is relevant. The spin dynamics then enters through the 
additional factor V±(t) d±V±(0) that describes the precession of a two-component spinor 
along the trajectory 7 ± in the given external fields E and B. In order to see this recall 
that V±(p,x), being an isometry from a subspace in C 4 to C 2 , projects four- component 
(Dirac) spinors to two-component spinors. These two components are the expansion coef- 
ficients of the projection to the eigenspace P^(p, a?)C 4 in the basis {ei(p, x), e^ip, x)} or 
{fi{p, x), /2(p, x)}, respectively The SU(2) matrix d± then propagates the two-component 
spinor along the trajectory ^ towards its end point. The propagated two-component 
spinor is then converted back into a four- component (Dirac) spinor by V±(t). One can thus 
view the combined dynamics of the spin and translational degrees of freedom as it enters 
the semiclassical propagator ( |4.1(J| ) to be of a mixed quantum-classical type: the transla- 
tional dynamics are classical and drive the quantum mechanical spin dynamics. However, 
it can be shown that expectation values of the spin operator, which is proportional to er, in 
the two-component spinors exhibit the well known Thomas precession of a 'classical' spin 
along 7 ± , see [24, 16, 25 1 for more details. Moreover, knowing the classical spin precession, 



one can recover the full spin dynamics, as described by the spin transport matrices d±, 
with the help of a certain dynamical and geometric phase fl6 ]. 



5 The trace formula with expectation values 

Having the semiclassical propagator (4.10) available, one could now proceed to derive a 



semiclassical trace formula, very much in analogy to Gutzwiller's original trace formula for 
non-relativistic, spinless particles [0, |2|. One only has to localise in energy first, in order 
to project out the essential spectrum of the Dirac-Hamiltonian. To this end one chooses 
a smooth function x £ Co°(^) that is supported in the interval (—mc 2 ,+mc 2 ) where the 
eigenvalues of H are located that one is interested in. Then x{H), defined by the functional 
calculus given by the spectral theorem, is a bounded and self-adjoint operator with pure 
point spectrum located in suppx; its eigenvalues are x(E n ). Hence the truncated time 
evolution operator U x {t) := e~~a tH x{H) has a Schwartz kernel with spectral representation 

K x ix,y,t) =J2x(E n )M*)^n(y)^ Ent (5.1) 

n 

in terms of eigenspinors ip n and eigenvalues E n G supp x °f H. Since K x (x,y,t) = 
x{H) K(x, y, t), where x{H) acts on the variable x, the leading semiclassical asymptotics 
of K x can immediately be concluded from ( j4.10|) . The only modification is provided by an 



additional factor of x(E^±), where E~± denotes the energy of the trajectory 7 . The result 



can then be used to derive a semiclassical trace formula; see ||T3|, |TIJ, where this procedure 
has been carried out. 



Here, however, our aim is to include the expectation values (ip n , Bip n ) of observables B 
as they have been introduced in section |3|. To this end one considers a suitably regularised 
trace of the operator B U x (t). The regularisation requires a test function p e C°°(1R) with 
compactly supported Fourier transform, p e C£°(IR), so that 

&x\fl '■=■*- I p{t)U x {t)dt (5.2) 

defines a bounded operator. The map p \— > Tr U x [p], where Tr(-) is the operator trace on TC, 
then yields a tempered distribution, if the sum ^2 n x(E n )p(E n /h) is absolutely convergent. 
The latter condition certainly holds when the eigenvalues do not accumulate in supp x- 
One hence should choose \ in such a way that possible accumulation points of eigenvalues 
lie outside of supp x, an d that x vanishes sufficiently fast towards accumulation points. 
Under this condition one can calculate the trace of the operator B U x [pe^ E ^], 

T * J p it)e f^BU x (t)dt = Y,x(E n )^ n ,Bi, n )p{^-^-) . (5.3) 



The right-hand side of ( |5.3|) already provides the spectral side of the trace formula we are 
aiming at. The semiclassical side of this trace formula can then be obtained from the left- 
hand side upon introducing the kernel K x for U x and calculating the trace in coordinate 
representation. 

Before one can perform this computation, however, one needs the Schwartz kernel for 
the operator B U x (t), or at least a semiclassical expression for it. The latter can be derived 
from the representation (|4.2|) of the time evolution kernel, when one employs the action of 
a Weyl operator D with symbol D(p, x) on functions of the form a(x) en^ x \ This reads 

(D a e*+) (x) = (D(V<f>(x), x) a(x) + 0(h)) e^ {x) , (5.4) 



see, e.g., [26]. In a first step one therefore has to apply the operator D = xifl) to K(x, y, t), 
which can be done with the help of ( |5.4j ) since according to the functional calculus of P2| 
x{H) is a Weyl operator with principal symbol D (p, x) = x{H(p, a?)). Then, in a second 
step, one repeats the procedure with the observable B. Thus 

= B (V^ix, y, t, £), x) xiH^V^ix, y, t, £), x)) (5.5) 



± 



«o 



x,y,t,^)ey ±{x ' y ' t ^ + 0(h) 



yields the leading semiclassical asymptotics for the integrand of the Schwartz kernel for the 
operator B U x (t) in an analogous representation to (|4.2|) . Calculating the trace on the left- 
hand side of (|5.3j ) in coordinate representation therefore requires, in leading semiclassical 
order, to evaluate the integrals 

/ / / ti(B (V x <f ) ± ,x)a±(x,x,t,£)) 
Jm.Jm.zJM.a (5.6) 

p(t)x{H ± {V x (t) ± ,x))e^ ± ^ x ^ )+Et) d^dxdt , 
10 



where tr(-) denotes a matrix trace. 

Since the integration cannot be performed in closed form, one employs the method 
of stationary phase. The connection of the phase with the generating functions S^ 1 then 
reveals that stationary points ($, st ,x st ,t st ) are such that (£ st ,x st ) lies on a periodic orbit 



7p with energy E and period t s t of the classical dynamics generated by H . One therefore 
concludes that up to terms of order h°° the expression (|5.3|) is given by a sum over the 
periodic orbits 7^ with periods T ± G suppp. We recall that the method of stationary 
phase also demands that at stationary points the phase be non-degenerate. Here this means 
that transversal to the (connected) manifolds of stationary points, which are at least one 
dimensional, the matrix of second derivatives of the phase with respect to (£, x, t) must 
have a (constant) rank of seven minus the dimension of the respective manifold. Let us 
explicitly mention two examples: 

1. The entire hyper-surfaces 

fi|:={(p,aj); H ± (p,x) = E} (5.7) 

of constant energy are manifolds of periodic points with trivial periods t s t = 0. Thus 
^e x {tst = 0} are five dimensional manifolds of stationary points. In this case the 
non- degeneracy condition means that E is a regular value for the Hamiltonians H^, 
i.e. (V p H ± (p,x),V x H ± (p,x)) ^ for all (p, x) G fl|. 

2. An isolated periodic orbit 7^ is non-degenerate, if its monodromy matrix M ±, i.e. 
the restriction of the linearised Poincare map to the directions transversal to the 
orbit, has no eigenvalue one. 

For simplicity, below we will always restrict attention to these cases, i.e. we require E to be 
a regular value for the Hamiltonians and assume that all periodic orbits with energy E are 
isolated and non-degenerate. The latter condition is, e.g., fulfilled by hyperbolic periodic 
orbits, where all eigenvalues of M + are different from one in modulus. However, we 
would like to stress that a semiclassical trace formula can be derived under more general 
conditions. One only needs a non-degeneracy condition for the manifolds of stationary 
points in order to be able to apply the method of stationary phase. 

For the explicit computation of the semiclassical side of the trace formula one should 
treat each connected manifold of stationary points separately. To this end we now introduce 
the partition of unity V • hj(t) = 1, with hj G C£°(R), under the integral (|5.6|) such that 
supp hj contains only one period T ' ± of a periodic orbit 7^ with energy E; in particular, 
h(T ±) = 1. Then each summand, labeled by j, corresponds to a contribution of exactly 
one (isolated and non-degenerate) periodic orbit. We, furthermore, reserve the label j = 
for the contribution of the manifolds f2^ for which t st = 0. The necessary computations 
are almost identical to the case of the trace formula without expectation values of an 



observable and thus we can direct the reader to [161 for the details. 



It turns out that the leading semiclassical order of the Weyl term, i.e. of the contribution 
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with j = 0, reads 
p(0) 1 



X(E) 



2tt (2nhy 



+ 



tr(P (p, x) P+(p, xj) 6(H + (p, x) - E) dpdx 
[ [ ti(B (p,x)Po(p,x))6{H-(p,x)-E)dpdx 



(5.* 



Introducing the notation 

d/i|(p, x) :-- 



1 



^ <*(#*(?, as) - E) dpdx 



vo\{Vt El 
for Liouville measure on Q, E , the leading order of the Weyl term can also be given as 

p(0) 1 



(5.9) 



X(E) 



2tt (2nhy 



vol(n+) / tr (P P +) d/4 + vol(^) / tr(BoP ") d^ 



. (5.10) 



Since the matrices Pq are projectors, the integrands in ( p,10| ) also read tr(P P { 



tr(P BqP ) =: trP , and as a short-hand we will also use trP for the integrals over 



In the case of an isolated and non-degenerate periodic orbit 7^ the evaluation of its 
contribution to the semiclassical side of the trace formula also proceeds in close analogy to 
the case without an observable, until it comes to the point where an integration over the 
primitive periodic orbit associated with 7^ is required. In the present situation the object 
that has to be integrated along the primitive orbit is the matrix trace of 

B {p, x) V±(p, x) d ± (p, x, T : Vl(p, x) 



'% 



(5.11) 



where the last three factors derive from the spin contribution to the semiclassical prop- 
agator ( f4.10| ). For the integration we parametrise the periodic orbit by (P ± (t), X ± (t)), 



< t < T ±, where T ± denotes the primitive period associated with 7^, and hence have 

ip ip " 

to evaluate ( |5.11D at (p,x) = (P (t),X (t)). The quantity that finally enters the trace 
formula then is a particular average of the principal symbol Po of the observable along the 
primitive periodic orbit, 



trP7 p 



1 



T* ± Jo 
ip 



■£ 



tr 



B (P ± (t),X ± (t))V ± (P ± (t),X ± (t)) 
d^PHtlX^T^vKP^X^t)) 



(5.12) 



dt . 



The integrand appearing in (|5.12j ) can be interpreted as follows: the right-most term 
projects a four- component (Dirac) spinor at (P ± (t), X ± (t)) to the local P ± -eigenspace 
P ( ^ = (P ± (t), JC ± (t))C 4 in a two-component representation. This two-spinor is then propa- 
gated once around the periodic orbit and subsequently converted back to a Dirac spinor 
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which is finally mapped by Bq. The integration then averages the result of this procedure 
over the primitive orbit related to 7^. 

Altogether, if all periodic orbits generated by the two Hamiltonians H ± with an energy 
E that is not a critical value for these Hamiltonians are isolated and non-degenerate, the 
semiclassical trace formula we have been aiming at reads 

' E„ — E s 



^2x(E n ) (if) n , Bij) n ) p 



h 



X(E) ^ j2^Ef ( V ° KQ+e) trB f + ™K^) tr V) +0(/T 1 ) (5 . 13) 
+ £x(£) ^ tr Z^ A lt e^ ±(S) (l + 0(h)) , 



2vr 

it 



where as usual the action of a periodic orbit is denoted as 

S^(E):= f pdx, (5.14) 



'7p 

and the amplitude associated with a periodic orbit, 

„# -if M ± 

i ± e 7 p 
4v± := - 777 , (5.15) 

7p |det(M + -l 4 )| 1/2 l * 

incorporates the monodromy matrix M ± and the Maslov index \x ± of that orbit. 

Let us remark that this trace formula reduces to the one derived in |I5*. IB"] in the case 



B = 1ft. Obviously, on the left-hand side of (|5.13|) the expectation value then disappears 



and in the Weyl term tr Bq = 2. In order to see what happens to the weights tr5 p 

of periodic orbits we notice that after a cyclic permutation under the trace in ( p-12| ) one 

can employ the fact that V±V± = 1 2 , and that trd±(p,x,T±) =: trd ± is independent of 

± 

(p,x); hence tr B Q p = trd ±. With these identifications the trace formula of fl5], |II| is 

hence recovered. 

As a final comment on the trace formula ( |5.13f ) let us mention that one can view the 

left-hand side of ( |5.13| ) as resulting from an application of the weighted spectral density 

d X:B (E) :=J2x(E n )(^n,B^ n )5(E-E n ) (5.16) 

n 

to the test function p (after a suitable shift of the variable). In this context the trace 
formula ( |5.13|) can be converted into a distributional identity for the weighted spectral 



density (|5.16|) , 

d x , B (E) = X(E) p^p (™l(Si+) tr Bf E + vol(fi;) tr B^) + 0{H' 2 ) 

+ *W h £ *** H ^' E) (1 + 0(K)) • <5 ' 17) 



2nh 

"1v 
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It is in this form that trace formulae are frequently presented (see, e.g., Gutzwiller's original 

work M). 



6 Semiclassical averages of expectation values 

The left-hand side of the trace formula ( |5.13| ) can be seen as a weighted superposition 
of the expectation values (ip n , Bip n ) corresponding to the eigenvalues E n in some interval 
about E of width proportional to H. This is so since p G C°°(R) is a test function with 
Fourier transform p G Cq°(M) C <S(1R). Thus p itself is, indeed, a Schwartz function, 
i.e. it is rapidly decreasing, and therefore essentially cuts off all terms corresponding to 
eigenvalues outside of some interval [E — huoi, E+fuv^]. Although the length of this interval 
shrinks to zero in the semiclassical limit, the number of eigenvalues it contains grows, since 
according to the Weyl term of the spectral density (compare ( |5.17| ) with B = In) the 



mean density of eigenvalues increases like h~ 3 . In the following we therefore want to study 
a semiclassical average of the expectation values (if; n , Bip n ) related to the Nj eigenvalues 
E n G I(E; h) = [E — friu, E + Hou], see ( |2.6| ). More precisely, we want to express the limit 



im n^T E ^Hn) (6.1) 



hm 
n 

' E n £l(E;h) 



in terms of classical quantities. 

It seems that the trace formula ( |5.13|) is the ideal tool that allows one to achieve the 



calculation of the limit QS.lp . For two reasons, however, the situation is not so straight 
forward: 



1. For the trace formula (|5.13| ) we assumed that all periodic orbits be isolated and non- 
degenerate. This is a condition which is, e.g., met by hyperbolic classical dynamics. 
However, hyperbolicity is a strong chaotic property which is not at all necessary for 
the result we are aiming at. In fact, a considerably weaker condition suffices: it is 
only required that the points on the energy shells £1^ that lie on periodic orbits are 
of Liouville measure zero. Any ergodic dynamical system shares this property, but 
even many integrable ones as well. 

2. For the test function p in the trace formula one cannot simply choose the characteris- 
tic function X[-u,u] of the interval [—a;, to], which would provide the sharp cut-off ( |6.1| ) 
in the summation over the eigenvalues, since X[-u,u] is n °t smooth and, moreover, its 
Fourier transform | smut is not compactly supported. It even decreases too slowly 
to make the sum over periodic orbits in (|5.13| ) convergent. 

Despite of these objections it is nevertheless possible to calculate the limit flS.lp in rather 
general situations. For the case of Pauli-Hamiltonians this has, indeed, been carried out 
in ||19| . However, to achieve this one does not employ the trace formula ( |5.13| ) directly, 
but rather starts off from (15. 31), which served as the basic relation from which the trace 
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formula was derived. Again one has to calculate the integral ( |5.6| ), but here we are only 
interested in the leading semiclassical order. By an inspection of the trace formula one 
realises that this is completely determined by the Weyl term, i.e. by the stationary points 
(£ st , x st , t st ) with t st = and (£ st ,x st ) G Q E . As a non-degeneracy condition we therefore 
only need that E is not a critical value for the two Hamiltonians H ± (p 1 x). For later 
purposes we, however, require this condition to hold for all E' G [E — e,E + e] with some 
e > 0. It is then ensured that the non-trivial time components t st ^ of stationary 
points cannot accumulate at zero (see, e.g., P^|). We stress that no assumptions on the 
periodic orbits have to be made apart from the condition that they form a measure-zero 
set on the energy shells. Under the integral ( |5.6| ) one now only needs a simple partition of 
unity, 1 = h (t) + (l — h (t)), that separates the stationary points ($, st ,x st ,0) from those 
with t st ^ 0. In order to estimate the latter contribution one exploits the fact that the 
points on periodic orbits, which are precisely the stationary points (£, st , x st ) with t st ^ 0, 
are assumed to be of Liouville measure zero on Cl E . As a consequence one obtains that 
the leading semiclassical order is solely determined by the Weyl contribution (see, e.g., 
|27|, |Tg], 0), so that 



(6.2) 

= X{E) 27 J^W ^° m+E) tT ^ E + V ° KQe) tr ^) + ° ( ^ 2) ■ 

In a next step one would like to get rid of the test function p and to replace it by a sharp 
cut-off. Under the conditions mentioned above, which led to the relation ( |6.2j ), this can 



indeed be done with the help of the Tauberian Lemma of p8| . The only point that still 
has to be clarified is that the Tauberian Lemma requires the expectation values of B to be 
non-negative. Since, however, B is bounded, this condition can be satisfied after shifting 



the observable by a suitable constant. In effect, the Tauberian Lemma of |28[ then allows to 
simply replace p by the characteristic function X[-u>,w]j an d its Fourier transform by | sin cut. 
We also recall that E was chosen to be somewhere inside the gap (—mc 2 , +mc 2 ) of the 
essential spectrum of H. Thus, for sufficiently small ft the interval I(E; K) is completely 
contained in (—mc 2 , +mc 2 ) so that the function x G C^°(M) that serves as to cut off the 
essential spectrum can be chosen such that x(E) = 1 = x(E n ) for all E n G I(E; h). Hence 

J2 <V>n, Hn) = ^ j^ (vol(ft+) tr B+ E + vol(fii) tr!^) + o(/T 2 ) . (6.3) 

E n £l(E;H) 71 { TV ) 

Furthermore, the choice B — 1% in ( |6.3|) yields a semiclassical representation for the 
number Nj of eigenvalues in the interval I(E; h), 

Ni = ^ 7^|p (vol(fi+) + vol(f^)) + o(h~ 2 ) , (6.4) 
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so that ( |6.3| ) and ( |6.4j ) finally allow to determine the limit ( |6.1| ), 

h -" Nl E n mE-M 2 vol(tij) + voi(na l ; 

As desired, the right-hand side of this relation involves only classical quantities. 

The result (6.5) can be converted into a statement about averages of the Wigner trans- 
forms 

W[ip n ](p,x):= [ e -* p - y i^(x-ly)®ipn(x+ly)dy (6.6) 

Jr 3 

of the eigenspinors ip n G I(E; h), since expectation values of Weyl operators can be ex- 
pressed as 

(ip n ,Bi; n ) = — i— / / tr(W{i} n ](p,x)B(p,x;h)) dpdx . (6.7) 

{inhy j r3 j R 3 

Here tr(-) denotes a matrix trace, whose presence prevents an immediate conclusion from 
( |6.5| ). To proceed nevertheless, we remark that if one replaces the observable B by its 
diagonal part Bd, see ( |3.4|) and (3.5), the right-hand side of ( |6.5|) remains unchanged. This 
immediately follows from the form ( |3.6| ) of the principal symbols of the diagonal part and 
the non-diagonal part, respectively. Thus, in particular, 



im 7T 52 i^n, B nd ip n ) = . (6.1 



hm 

h 

E n &I{E;h) 



We can therefore now restrict attention to observables B that are Weyl quantisations of 
symbols B = B 0i d = B^ + B . Starting with a quantisation of, e.g., Bq one splits this 
matrix valued symbol into a sum of contributions where all matrix entries but one vanish, 
so that now one essentially deals with scalar operators and symbols. Applying ( |6.5| ) to 
these then finally yields 

E n eI(E;h) ,g g s 

_ 1 yol(g|) 6(H+(p, x) - E) P+(p, x) + vol(^) S(H~(p, x) - E) P ~(p, x) 
2 vol(fi+) + vol(^) 

Here the convergence as ft —>■ has to be understood in the sense of distributions, i.e. 
after an application to (matrix valued) test functions. In (|6.5| ) the latter are given by the 
symbols of Weyl operators. 

Let us add a few remarks on the result ( |6.5|) . Relations of this kind, which express 
semiclassical averages of quantum observables in terms of classical averages are known 
as (weak versions of) Szego limit formulae p9[ . They are well known for Schrodinger- 
Hamiltonians (see, e.g., |2"7fl ) and have recently been established for Pauli-Hamiltonians 
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|H| . In the present case the classical side contains the two contributions that arise from 
projections of the classical observable B to the two eigenspaces of the symbol H and 
hence correspond to 'particles' and 'anti-particles', respectively. The relative weights of 
the two terms are fixed by the relative volumes of the corresponding energy shells. It often 
occurs that for a given energy E only one of the energy shells is non-empty, and then only 
one contribution appears on the right-hand side of ( |6.5| ). The only exceptional situations 
are those where the Klein paradox occurs, i.e. where a tunneling between particle and 
anti-particle states is possible. The corresponding mixture is then accounted for by the 
presence of both terms on the classical side of ( |6.5|) . Due to the relation (|6.8| ), however, 
semiclassically transitions from particle to anti-particle states, or vice versa, play no role. 
On the right-hand side of ( |6.9|) this effect is caused by the fact that only terms proportional 
to the projectors P^ occur so that the semiclassically averaged Wigner transforms are 
block-diagonal with respect to the particle and anti-particle subspaces. 

Obviously, Szego limit formulae of the kind Q6.5| ) contain considerably less information 
on expectation values than trace formulae like (|5.13| ). They demonstrate that the semi- 
classically leading term of the weighted spectral density ( |5.17| ) is the semiclassical average 
of the observable multiplied by the leading Weyl-term of the spectral density. Fluctuations 
of expectation values about the semiclassical average are then described by the sum over 



periodic orbits in ( |5.17|) and therefore depend on classical properties more sensibly than 



the average itself. However, the trace formula (|5.13|) does not allow to draw conclusions 
about individual expectation values since on its spectral side it contains the test function p, 
whereas on its semiclassical side the Fourier transform p appears. Thus, by Fourier duality, 
shrinking the (effective) support of p results in a growing (effective) support of p, finally 
leading to a divergent, and thus uncontrollable, sum over periodic orbits. The smallest 
support of p that can, indeed, be handled leads to the semiclassical average flfc>.5|). How- 
ever, in the same way as the Gutzwiller trace formula is not suited for the representation 
of individual eigenvalues, but enables a semiclassical analysis of spectral correlations (see, 
e.g., 0, 0), the trace formula ( |5.13| ) can be used to determine correlations of expectation 



values as, e.g., in |31 | 
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